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Abstract 

Abstract: We show that the kinetic approach to statistical mechanics per¬ 
mits an elegant and efficient treatment of fractional exclusion statistics. By 
using the exclusion-inclusion principle recently proposed [Phys. Rev. E 49 , 
5103 (1994)] as a generalization of the Pauli exclusion principle, which is 
based on a proper dehnition of the transition probability between two states, 
we derive a variety of different statistical distributions interpolating between 
bosons and fermions. The Haldane exclusion principle and the Haldane-Wu 
fractional exclusion statistics are obtained in a natural way as particular cases. 
The thermodynamic properties of the statistical systems obeying the gener¬ 
alized exclusion-inclusion principle are discussed. 
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Following Haldane’s formnlation of a generalized Panli exclusion principle, many 
papers have been recently devoted to the study of fractional exclusion statistics by interpo¬ 
lation of bosonic and fermionic distributions . There is an intrinsic connection between 
these fractional statistics and the interpretation of the fractional quantum Halleffect [Q and 
anyonic physics Murthy and Shankar generalized the Haldane statistics to infinite 
dimensions and showed 0] that the one dimensional bosons interacting through the two- 
body inverse square potential Vij = g{g — l){xi — Xj)~‘^ of the Calogero-Sutherland model 
obeys fractional exclusion statistics in the sense of the Haldane interpretation. Using the 
Thomas-Fermi method, Sen and Bhaduri [Q considered the particle exclusion statistics for 
the Calogero-Sutherland model in the presence of an external, arbitrary one-body poten¬ 
tial. Isakov 1^ extended the Calogero-Sutherland model in the case of particles interacting 
through a generic two-body potential Vtj and derived the related exclusion statistics. Ra- 
jagopal obtained the von Neumann entropy per state of the Haldane exclusion statistics. 


Nayak and Wilczek studied the low-temperature properties, the fluctuations and the 
duality of exclusions statistics. Murthy and Shankar ^ also computed, for quasi-particles, 
in the Luttinger model, the exclusion statistics parameter g related to the exchange statis¬ 
tical parameter a of the quantum phase e*™, and showed []^ that the parameter g can be 
completely determined by the second virial coefficient in the high temperature limit. 

Haldane dehned the statistics g oi a. particle by 


9 = 


^JV+AJV 

AN 


( 1 ) 


where is the single-particle Hilbert-space dimension, when the coordinates of — 1 
particles are kept fixed. Thus the dimension of the Hilbert space for the single particle 
states is a hnite and extensive quantity that depends on the number of particles in the 
system. One can choose this dimension as 


d^ = K-g{N-l) 


and the statistical weight or degeneracy factor is 


( 2 ) 
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w = 


(d.+iV-l)! 


( 3 ) 


7V!(d,-l)! • 

This equation is a simple interpolation between the number of ways of placing N identical 
bosons or fermions in K single-particle independent states confined to a finite region of 
matter. The expression of the mean occupation number n = N/K has been obtained in an 
implicit form by Wu 0]. 
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n = 


w{E) + g ’ 


( 4 ) 


with w satisfying the relation 


w{Ey [1 + w{E)Y-3 = . 


( 5 ) 


In the special cases g = 0 and g = I, Eq-(4) yields, respectively, the Bose-Einstein (BE) and 
Fermi-Dirac (ED) distributions. 


Ilinski and Gunn ||^ criticize the identity of relations (1) and (3) of Haldane and, 
using (2), derive a different W constructing a statistical mechanics which is not in complete 
agreement with Wu statistical mechanics [Q. 

Acharya and Swamy |]^ and Polychronakos have studied, in addition to the Haldane 
exclusion statistics, new fractional statistics with appealing features as positive probabilities 
and analytical expressions for all termodynamic quantities. We wish to recall the first work 
on intermediate statistics published by Gentile at the beginning of 1940 [|T^ . 

In a previous work [|^ we have considered the kinetics of particles in a phase space of 
arbitrary dimensions D, obeying an exclusion-inclusion principle. We obtained the statistical 
distributions of the particle system as the stationary state of a non-linear kinetic equation. 
A crucial point of this formalism is the definition of the transition probability which can be 
written in various forms, containing the effects of the inclusion-exclusion principle through 
the distribution function n = n{t,v), which is the mean occupation number of the state v. 

Setting 7r(t,v — u) the transition probability from the state v to the state u, the 
evolution equation of the distribution function n{t, v) can be written as 
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dn{t, v) 
dt 


J [7r(t, u ^ v) — 7r(t, V —> u)] 


u 


( 6 ) 


In ref. |16] we postulated the following expression of the transition probability : 


7r(f, V u) = r(f, V, V - u) ip[n{t, v)] u)] , 


(7) 


where — u) is the transition rate, (p[n(f, v)] is a function depending on the occu¬ 

pational distribution at the initial state v and 'ilj[n{t,u)] depends on the arrival state. The 
function (p{n) must obey the condition (p(0) = 0 because the transition probability is equal 
to zero if the initial state is empty. Furthermore, the function 'ip[n) must obey the condition 
"0(0) = 1 because, if the arrival state is empty, the transition probability is not modihed. If 
we choose ^{n) = n and 'ip{n) = 1 we obtain the standard linear kinetics. The function 'ip{n) 
defines implicitely the inclusion-exclusion principle enhancing or inhibiting the transition 
probability. 

In this letter we use the kinetic approach outlined above to generate new fractional 
exclusion statistics. The approach is appropriate both for non-interacting and interacting 
particles. In the following, we derive general expressions linking the main statistical and 
thermodynamic quantities to the transition probabilities of the kinetic theory. We obtain 
implicitely a fractional statistics which interpolates between BE and FD distributions in a 
single-particle Hilbert space whose dimension is an arbitrary function of n. We show that the 
Haldane-Wu statistics can be derived from this one by considering Brownian particles and 
by demanding that the Hilbert space dimension be a linear function of n. Another particular 
case within the family of statistics here introduced yields a statistics interpolating among 
BE, FD and Maxwell-Boltzmann (MB) distributions. 

Let us consider Eqs. (6) and (7) and limit the discussion to the one-dimensional velocity 
space (the extension to a D-dimensional phase space is straightforward). We limit ourselves 
to first neighbors interactions, this being equivalent to truncate the Moyal expansion given 
by Eq.(7) of Ref. |]^ at the second order and we obtain the following, generalized, non-linear 
Fokker-Planck equation 
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(8) 


dn{t,v) d 

dt dv 


J{t,v) + 


dD{t,v)' 


ip{n)'ilj{n) 


+D{t,v) 


dv 

d'ip{n)\ dn{t,v) 


dv 


dn dn 

J{t,v) and D{t,v) are the drift and diffnsion coefficients, respectively, and are given by the 
hrst and the second order moments of the transition rate []^. Eq.( 8 ) is a continnity eqnation 
for the distribntion fnnction n = n{t,v) 

dn(t,v) dj(t,v,n 


dt dv 

where the particle cnrrent j = j{t, v, n) is given by 


= 0 


(9) 


j = -D 


ip{n)'d{n) 

{Any 


de 


dn 


IT. 


dv 


dv 


( 10 ) 


We have defined 


d 


{Any = <^ — log 


y{n) 


n ^ -1 


( 11 ) 


and 


dv 


D{t,v) 


J{t,v) + 


dD{t, v) 
dv 


( 12 ) 


The fnnction e = e{v) is an adimensional single particle energy defined np to an additive, 
arbitrary constant. This energy is appropriate both for non-interacting and interacting 
particles. Nayak and Wilczek |]^ have examined the particular case e{v) cx v^ with I integer. 
The case I = 2 corresponds to Brownian particles, the drift and diffusion coefficients being 
given by 


J = B = ^ . (13) 

The quantity 7 is a dimensional constant and e is given by: e = P{E — fi) with E = 
being the kinetic energy. (3 = 1/T {kg = 1) is the inverse of the temperature and /i is the 
chemical potential. 

In stationary conditions, the particle current vanishes: j{oo,v,n) = 0. Eq.(lO) becomes 
a omogeneous first order differential equation, which can be integrated easily: 
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'ipin) 

ip{n) 


exp(e) 


(14) 


When the functions (p{n) and 'ip{n) are hxed, Eq.(14) gives the statistical distribution n = 
n{e) of the system. 

We stress that, if we use as variable the single particle energy e, it is easy to verify that 


the above-dehned function (An)^ =< > — < n >^, is equal to the second order density 


fluctuation 17 


(An)2 = -— . (15) 

This relation is very important because it reveals that the mean fluctuation of the occupation 
number is the crucial quantity that determines the equilibrium distribution. 

The statistical distribution n(e) can be obtained by the maximum entropy principle, 
hxing the total number of particles and energy of the system and using the standard Lagrange 
multiplier method. Setting S{n) = S{N)/K the thermodynamic limit of the entropy per 
state, we have 


L[S(n,)-eii| = 0 , 


( 16 ) 


From Eqs. (14) and (16) it is possible to obtain the entropy in terms of the functions 
(p{n) and '^(n) as 


tfjin) 

ip{n) 


= exp 


dS{n) 

dn 


(17) 


or in terms of (An) 


(An)2 = — 


d^S 

dm? 


-1 


( 18 ) 


By using the Boltzmann principle with the identification of the entropy S = logW, it is 
possible to find the relation between W and the functions (^(n) and 'ip{n) in the following 
expression 


■^(n) 

(p{n) 


1 dW{n) 
>V(n) dm 


( 19 ) 
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It is well known that the partition function Z per state 



( 20 ) 


is related to the mean density n by means of 

n = -^logZ . (21) 

Taking into account Eqs. (11), (15), (20) and (21), we may write the thermodynamic 
potential hi as a function of 93 (n) and 'ip{n) 

pn= f n-^log dn , ( 22 ) 

J on ^{n) 

and the density fluctuation as 

(23) 

Now let us consider the case (p{n) = n; for the function 'ip{n) we choose the particular 
form 

ip = d sn(^n + dy~^sn . (24) 


{Any = —n 




where the function d = d{n) must satisfy the condition (i(0) = 1. In this case, Eq.(14) which 
dehnes the statistics, becomes: 


dd 

dn 


[n + d) 


1+1 


= ne 


The function W becomes 


and the entropy 


>V 


{n + dy~^'^ 
n" d‘^ 


(25) 


(26) 


S = —n[w log w — (1 + w) log(l + tc)] , 


(27) 


where w = d/n, so we hnd the general form of the von Neumann entropy per state |TO|. The 
thermodynamic potential hi becomes 
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( 28 ) 



(An)^ = 


nd{n + d) 

, dd\^ , , ,, d^d , n + d 

d-n—\ - nd[n +d)\og —— 
on / on^ d 


The quantity d is an arbitrary function and we can use it to define a family of exclusion 
statistics. If we select a particular form of d, we observe that this one can depend not 
only on n but also on a parameter g: d = d{g,n). In this way Eq.(25) dehnes, varying 
g E [0,1], a fractional statistics. If we require that the statistics interpolates between BE 
and ED distributions we must set the two conditions: d{0, n) = 1 and d{l, n) = 1—n. These 
two conditions imply that the single-particle Hilbert-space dimension be that of the Bose 
space when g = 0 and that of the Fermi one when g = 1. Let us call S{N) = KS{n) and 
W{N) the entropy and statistical weight of the system of N particles lying in K states. 
From the Boltzmann principle S = loghE we obtain W = It is easy to see that 

W, given by Eq.(26), is the thermodynamic limit of the statistical weight given by Eq. 
(3), where, this time, = d{g,N,K) is the single-particle Hilbert- space dimension, is 
an arbitrary function and admits as thermodynamic limit d{g,n) = d{g,N,K)/K. If we 
consider a statistics interpolating between bosons and fermions, the conditions to which 
d = d{g, n) must satisfy when = 0,1 can be derived in the thermodynamic limit from the 
two conditions d(0, N,K) = K and d(l, N, K) = K — N + 1. 

The Haldane-Wu choice of d given by Eq.(2) implies 


d = 1 — g n 


(31) 


and requires that ip{n) be given by 



^(n) = |1 - jn]’[1 + (1 - t/)n]‘ ® 


(32) 


Equation (25) becomes 


[1 — gn\^ [1 + [1 — g)nY ^ = ne 


(33) 


and we obtain, in the case of Brownian particles, the statistics introduced by Haldane 
and by Wu [Q. The partition function Z and the density fluctuation are given by 

l + (l-^)n 


Z = 


1- gn 


(34) 


(An)2 = n (1 — gn) [1 + (1 — g)n] . 

As a second example of fractional statistics we consider the statistics dehned by: 

1 


(36) 


n = 


(36) 


studied extensively in ref. 


exp(e) — K ’ 

For K = — 1, 0 and 1 one obtains the FD, MB and 


BE statistical distributions, respectively. In ref. [^] it is shown that the parameter k is a 
function of the exchange statistical parameter a appearing in the quantum phase e*™. This 
statistics is generated from 


'ip{n) = 1 + K,n . 


In this case the partition function can be written as 


(37) 


Z = (1 + KIl)'''- , 


(38) 


and it is easy to verify that the density fluctuation (An)^ and the entropy S are given by 


(An)2 = n(l + Kn) , 


(39) 


5 = — (1 + Ku) log(l + Kn) — n logn , 

K 


(40) 


while the statistical weight W becomes 
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w = 


(1 + 


n" 


(41) 


and can be obtained as the thermodynamic limit of the statistical weight; 


W =\k 


N 


K 

\k\ 




W = 


K 


N 


N\ 


(— + N-l]\ 
W = K^^ 


1 ! 


. fK 
N\ — 
V K 


K < 0 , 


K, = 0 , 


K > 0 . 


(42) 

(43) 

(44) 


In conclusion, we have shown that the effects of the exclusion-inclusion principle can be 
taken into account in the particle kinetics by means of a proper dehnition of the transi¬ 
tion probability. We have derived a non-linear evolution equation which admits as steady 
solutions the exclusion statistical distributions. We have linked the main thermodynamic 
properties of the system with the function ■^(n) which generates the exclusion statistics. As 
hrst application of the theory we have considered a family of exclusion statistics in which 
the single-particle Hilbert-space dimension is an arbitrary function of the mean occupational 
number and contains, as particular case, the Haldane-Wu statistics. As second application 
we have considered an exclusion statistics interpolating among FD, MB and BE distribu¬ 
tions. 
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